Triply-periodic Smectics 
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Twist-grain-boundary phases in smectics are the geometrical analogs of the Abrikosov flux lattice 
in superconductors. At large twist angles, the nonlinear elasticity is important in evaluating their 
energetics. We analytically construct the height function of a tt/2 twist-grain-boundary phase in 
smectic-A liquid crystals, known as Schnerk's first surface. This construction, utilizing elliptic 
functions, allows us to compute the energy of the structure analytically. By identifying a set of 
heretofore unknown defects along the pitch axis of the structure, we study the necessary topological 
structure of grain boundaries at other angles, concluding that there exist a set of privileged angles 
and that the n/2 and n/3 grain boundary structures are particularly simple. 

PACS numbers: 61.30.Jf, 02.40.-k, 61.72.Mm, 61.72.Bb, ll.lO.Lm 



I. INTRODUCTION 

Because of their stability and quantization, it is natu- 
ral to regard topological defects as independent degrees 
of freedom in systems with broken symmetry. Perhaps 
the most storied example is the Abrikosov phase of type 
II superconductors 0. There, in the London limit, vari- 
ations in the phase of the macroscopic wave function can 
be decomposed into a "smooth" part and a singular part, 
the former being analogous to spin waves. The singular 
component represents the topological defects - those lines 
(in the three dimensional case) around which the phase 
slips by 27r. The spin waves can be removed from the the- 
ory resulting in an effective theory of repulsive vortices 
. Being the ideal proving ground for the study of bro- 
ken symmetries, liquid crystalline states provide us with 
a variety of line-like and point-like defects which, because 
of the anisotropy of the surrounding medium, can enjoy 
orientationally dependent interactions. With the flux- 
lattice in mind, we consider the smectic-A liquid crystal 
which is a close analog of the superconductor ■ Indeed, 
both the superconductor and the smectic-A phase have 
a complex scalar order parameter, representing macro- 
scopic phase ordering of the Cooper pair wavefunctions 
and the one-dimensional periodic density modulation, re- 
spectively. The superconducting order parameter is min- 
imally coupled to a gauge field, while the smectic-A order 
parameter is coupled to the director modes of a nematic 
phase which exists at higher temperatures: 



F = -FsmA[V', n] + Fprankin] 



(1) 



where 



(2) 
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and -Fprank IS the standard Frank free energy for a (pos- 
sibly chiral) nematic: 

i^Fra„k[n] = J rf^x {Xi (V • n)' + [n • (V X n) + kof 
+K3 [n X (V X n)]' } (3) 

and where n is the unit director field. Though deceptively 
similar to, for instance, a gauge-fixed Landau-Ginzburg 
theory the differences are profound. Fluctuation effects 
are famously more complicated at both the transition 
and in the ordered state 0, 0, 0] ■ 

The Landau-Ginzburg-De Gennes theory predicts a ne- 
matic to smectic-A transition when r < 0. In the or- 
dered state I (■0) I 7^ and the gradient term in (0) favors 
states for which V?/' — iqontp. Because n is a unit vec- 
tor, the phase of ip necessarily varies in space. In the 
ground state, if we take the nematic to order along z, 
then tp = \ip\e'^'^°^ . In the "London" limit, we consider 
fluctuations only in the phase of ip and not \ip\. Writing 
ip = {ipl ex.p{iqQ[z — u{x, y, z)]}, fi = z — Sn, and expand- 
ing FsaiA to quadratic order in the fluctuations u and Sn, 
we have: 

F ~ \j dhl^Bid^uf + B{\I i_u- 5nf 

+ Ki (V_L • 6nf + K2 (V_L X 5nf 
+ K3 {djnf + 2K2koWi_ X dn I (4) 



where 5n is the projection of the director in the xy- 
plane (appropriate for small director fluctuations), B = 
2(7q|i/'P, and ko is a pseudoscalar which would set the 
cholesteric pitch. Even in this quadratic theory, we can 
see that the transverse mode of 6n decouples from the 
(Eulerian) layer displacement u{x,y,z). The gauge-like 
coupling of V^M to Sft sets Sn = V_lm at long distances. 
As a result, the transverse modes of the director are at- 
tenuated at a length scale A = y^K/B where K is on 
the order of the Frank elastic constants. Known as the 
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Meissner effect in superconductors this result shows that 
smectic order excludes twist, xSn. 

Departing from the quadratic approximations of 
we can make this observation more general. In the full 
nonlinear theory, ( IV'I ) is still nonvanishing when r < 
0. In the London limit we write ip = \tp\e^'^°^ where 
$(a;, y, z) is a phase field and the mass density is p(x) — 
Pq + ji/)! cos[(7o*i'(x)]. It follows that the smectic layers 
sit at the density peaks defined by qo^{x) — 27m for 
n E "Z, i.e. level sets of $(x). Again, the gradient term 
in requires that V<I> — n for constant so at long 
distances the director is perpendicular to the level sets 
or, in other words is parallel to the unit layer normal 
N = VfE'/lVf&l. However, the twist order parameter is 
k = n • (Vxn), and it is straightforward to check that 
if n oc V$ then k = 0. Thus, the smectic layers are 
incompatible with twist^. 

In order to relieve the incompatibility between smectic 
order and twist, it is necessary to allow \ip\ to vary and, in 
fact, vanish at isolated points to form topological defects. 
Indeed, the competition between molecular chirality and 
the existence of layers leads to the celebrated twist-grain- 
boundary (TGB) phase (olITS. This phase is the analog 
of the Abrikosov flux lattice P|, with screw dislocations 
replacing flux lines and molecular chirality replacing the 
magnetic field. In smectics-A, however, there is an addi- 
tional complication: the coupling of geometry to elastic- 
ity requires that the screw dislocations and layers rotate 
together. For small angles of rotation, the layer structure 
can be approximated using linear elasticity Q J., ,11). 
The underlying rotational invariance of the compression 
strain, however, necessitates certain essential nonlineari- 
ties which have profound effects on t he g round state en- 
ergetics and layer displacements [T^ Il3l ll^. For large- 
angle twist-grain boundaries, seen for instance in bent- 
core systems [T5l| . the nonlinearities become important 
and we are forced to confront the full nonlinear elastic- 
ity. Because of the difficulty in systematically removing 
the unit director n from the theory, we go immediately 
to a rotationally-invariant free energy in terms of the 
phase field ^(x). Again, the gradient coupling in (0) sets 
V<I> = n, so we know that the compression strain van- 
ishes when IV^I = 1 or, equivalently, when N • V^" = 1, 
which indicates that the layers are spaced by one period 
along the layer normal - i.e. equally spaced layers. We 
write the compression strain in terms of both <I> and the 
Eulerian displacement field m(x) = z — $(x): 



1- (v$)-^ 



dzu - ^ {\7uf 



(5) 



The factor of ^ is introduced so that in the linearized 
strain, u^z — dzU, the standard result. Note that the 
nonlinear term in ((Sj is required by rotational invariance 



and is not merely an anharmonic correction to the elas- 
ticity. It is responsible for strong and subtle corrections 
to linear elasticity |0, |3 ■ The bending energy is 
inherited from the Frank free energy and is simply V • N. 
We note, however, that this is precisely twice the mean 
curvature of the layers H = • N. Finally, we 
write the full nonlinear free energy as 



F 



B 
1 



(V<5)' 



(6) 



where B is the compression modulus and = Ki/B is 
the "splay penetration length" . The challenge is to find 
ground states which minimize © with the appropriate 
boundary conditions. At infinity we take the layers to 
be perpendicular to the z-axis so that limx-»oo — 
z. There are other boundaries, however, namely lines 
corresponding to topological defects where $ changes by 
27r/go- 

Recently, we have developed an approach to study 
TGB phases with 7r/2 grain boundaries 0|. Our con- 
struction relies on a duality: 7r/2 twist-grain boundaries 
can be equivalently constructed from a sum of disloca- 
tions with burgers vector 6, or a sum of dislocations ro- 
tated by ninety degrees with burgers vector —h. Using 
this duality, a closed-form expression for the height of the 
smectic layers ensues and generates Schnerk's first sur- 
face (see Fig. nj, one of a class of mathematical surfaces 
formed by summing individual screw dislocations. In this 
paper, we elaborate on our construction and the energet- 
ics of Schnerk's first surface. We then identify a third set 
of defects that lie along the pitch axis. We use this set 
of previously unknown defects to study rotation angles 
smaller than 7r/2 and identify a privileged set of angles 
with particularly simple structures and correspondingly 
low energies. 

The building blocks of any TGB phase are the screw 
dislocations - when properly put together these effect a 
twist on the smectic layers without requiring the com- 
plete disappearance of the smectic order along a two- 
dimensional wall. We write the screw dislocation in terms 
of the phase field as: 



$ft(x) = z - — tan-M ^ 
2tt \x 



(7) 



^ The converse holds as well. Namely, if n ■ (Vxn 
surfaces can be found with n as their unit normals 



then 



It is easy to verify that this phase field is an extremum of 
equation ©. In addition, for b = ±2, this surface is iden- 
tical to a relic minimal surface [H = 0) jl^, the helicoid, 
and is harmonic in both three dimensions (V^$/i = 0) 
and in the xy-plane (V^$ = 0). Minimal surfaces have 
long been used as ansatzen for smectic structures be- 
cause they are global minima of the bending energy |l9j | . 
Nonetheless, the interplay between the compression and 
bending energies ensures that true minimizers of equation 
©, with more complicated topologies than the hehcoid, 
will neither have a harmonic phase field nor be minimal. 

In this paper, we will take the opposite tack and con- 
sider, not minimal surfaces, but harmonic $ formed by 
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dislocation of Burgers scalar b is given by 



FIG. 1: Schnerk's first surface, with charge +2 and —2 screw 
dislocations. Note that the —2 dislocations lie at the center 
of the rectangle made by the adjacent +2 dislocations. We 
choose e = Tp andk'^ ^ -0.03033 so that K'{k)/K{k) =2-i. 



taking arbitrary sums of parallel screw dislocations. If wc 
take the dislocations to have their defect axes along the 
z direction with positions = {xi,yi) in the sy— plane, 
then we write 



$(x) =72:-^$,,(x 



(8) 



As we will demonstrate, these sums will allow us to gen- 
erate models for tt/2 TGB phases. Furthermore, there 
is an unexplored but close connection between surfaces 
generated by equation © and certain triply- and doubly- 
periodic minimal surfaces. 

This paper is organized as follows. In sectional we 
review the nonlinear theory of large-angle twist-grain 
boundaries, and establish our notation. In section IIIII 
we present our theory for the 7r/2 TGB structure based 
on Schnerk's first surface. In section Hvl we discuss the 
properties of Schnerk's first surface. In scction[V] we dis- 
cuss the possibility of a moire phase of screw dislo- 
cations in TGB structures of arbitrary angle and discuss 
a potential lock- in mechanism for certain large angles. 
Finally, in section IVIl we summarize our results. 



II. TWIST-GRAIN BOUNDARIES AND THEIR 
DUALS 

Here we review the construction of a single twist- 
grain boundary with arbitrarily large angles. To simplify 
our analysis, we assume the smectic layers have normal 
N = z in the absence of dislocations and introduce the 
complex variable w = x + iy for the coordinates along 
the plane normal to N. In this notation, a single screw 



— Imln(u;). 
27r 



(9) 



We introduce the constant 7 here, which is set by enforc- 
ing the boundary condition that the compression strain 
vanish at infinity. In this simple case of one screw dislo- 
cation, we set IV^I — > 1 at infinity, implying that 
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1 



lim 



2tt\w\ 



1. 



(10) 



As we shall see, when considering doubly- and triply- 
periodic smectics, this boundary condition will depend 
on the orientation and position of the screw dislocations. 

A single twist-grain boundary can be decomposed into 
a set of dislocations with spacing id and positions w+idn 
for integers n Though we know of no minimizer 

of equation © with this property, progress can be made 
by considering the sum 



<I>row(2:, w) = 7Z - lmln{w + £dn 

Zn 



n— — 00 



(11) 



Utilizing the infinite product sin w 
follows that [21II 



wUn=^o (1 



b f TTW 

3'row = 7^ - — Imlnsm — 
In \ Id 



(12) 



To set 7 now, we must go out along the y-axis only, 
since the x-axis contains the grain boundary. Setting 
|V$row| — > 1 as 2/ ^ ±00, we have 



b 
2fd 



(13) 



The layers are flat at y = — 00 with layer normal given 
by N_ = ^2?7^ ^'^'^ y ~ layer normal is 

N+ = + 72. This gives an overall angle of rotation 

a^2s\n-^[b/(2ld)] Hill. 

This structure has the same topology as Scherk's first 
surface |l8j . a minimal surface which has been used 
as a coniectured structure for TGBs in diblock copoly- 
mers |l3- Surprisingly, this sum and Scherk's first sur- 
face are merely uniform dilations of each other: the level 
sets of $mw(a;, y cos(a/2), z) are identical to Scherk's first 
surface when b — 2. Furthermore, the rescaled sur- 
face is minimal for any b despite the fact that the mean 
curvature is nonlinear. It is reasonable to conjecture that 
the true smectic layer geometry falls somewhere in be- 
tween the exact sum of dislocations and Scherk's surface, 
and this expectation has been borne out by recent nu- 
merical studies of both lyotropic g2j and thermotropic 
smectics |23|. Numerically, $row (without rescaling) is a 
slightly better than Scherk's surface as a model for the 
layers. An important caveat, however, is that at large an- 
gles the director decouples from the layer normal near the 
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dislocation cores |23 , and it would be necessary to rein- 
troduce director modes to make contact with simulation. 
Nevertheless, approximate analytical models are useful 
guides for our understanding and, in this case, yield in- 
sights into the TGB smcctic structure. 

The coordinates {x,y, z) of the level set 'I'low — sat- 
isfy M 

tan(27r7z/6) tan(7ra:/^<i) = tanh^Try/id). (14) 

Equation H14|l possesses a hidden symmetry: if we rotate 
the entire structure by 7r/2, interchange {x, z) (z, —x), 
and simultaneously take 6 — s- —6, the equation is invari- 
ant. For a general rotation angle a, this transformation 
exchanges and 6/(27) in equation (|14|l and changes 
the rotation angle from a to tt — a. This allows us to 
view the level sets of a single TGB of any angle a as be- 
ing constructed either from parallel defects with burgers 
scalar h (along z) or from parallel defects with burgers 
scalar —b rotated by 7r/2 (along x). 
On the level sets of $row, we calculate 

27r7 cos^(27r7z/6) J 27r7tan(7rx/^d)z 
row = - . l{-Kx/td) \ b cos^{2n jz/b) 



+ 



b tan( 

TT tan(27r7z/5)a; 
idcos^inx/id) ^ 4 cosh^(7ry/4) 



Try 



(15) 



For a = 7r/2, when l-n^jh — i:/(,d, it is clear from equa- 
tion (|15|) that {x, z) (z, —x) preserves the unit normal 
vector N. The magnitude of V$row, however, is not pre- 
served by this transformation. To see why this must be 
the case, notice that |V$row| diverges at the core of de- 
fects along the z axis, even though we have now identified 
an additional set of helicoids along the x axis. Though 
rotations by ninety degrees preserves the level sets of the 
surfaces, as indeed they must, they need not preserve the 
structure of the phase field between those level sets and, 
in particular, do not preserve the orientation of the defect 
cores. We conclude that only one set of helicoids need 
have cores. However, the total compression energy and 
bending energy remains unchanged by a ninety degree 
rotation - the direction of the defect cores does not alter 
the energy. 



III. SCHNERK'S FIRST SURFACE 

A. The Tv/2 TGB phase structure 

A 7r/2 TGB phase consists of twist-grain boundaries 
with defects alternating along the z and x directions 
(we assume y is the pitch axis) and with rotation an- 
gle a = 7r/2. Finding an analytical expression for such 
a structure from which the energy can be calculated is a 
daunting task 12]. However, the degeneracy in how we 
identify the defects in a single TGB structure allows a 
significant simplification in the structure of this phase. 



By applying the appropriate ninety degree rotations, all 
of the dislocations can be rotated to be parallel to the 
z axis, leading to a structure of alternating TGBs such 
as the one shown in Fig. ^ This allows us to compute 
analytically the level sets of a tt/2 TGB phase by sum- 
ming parallel, alternating TGBs separated by a distance 
£b- This sum is given formally by 



tgb(x) =7Z-— Im ^ (-l)™lnsin 



TTW 

47 



(16) 

where r is a complex number that generates the ap- 
propriate translation between grain boundaries. To en- 
sure that each grain boundary rotates the layers by 7r/2, 
we might set, using the results of the previous section, 
b/{2£d) = sin(7r/4) = 1/V2. However, this relation re- 
quires a modification owing to the collaborative effect of 
the adjacent grain boundaries. We put off the details of 
setting the rotation angle to later. For now, we keep £d 
and b as free parameters. 

We may rewrite the infinite sum in $tgb Hlfijl as 



^'tgb(x) = 7z - — Im Ine(w) (17) 

ZTT 



where 



&iw)^ n 



sin {iTw/id + nmr /2) 
sin [-Kw/id + {m+ l)7rr/2] 



Because $tgb is a sum of harmonic functions in x and 
y, it follows that lnO(ti;) is analytic as is /(w) = e®'^™-'. 
Moreover, f[w) is doubly periodic and it thus follows that 
f{w) can be represented in terms of elliptic functions. In- 
deed, through an appropriate rescaling of x and y, f{w) 
shares all the poles and zeroes of the Jacobi elliptic func- 
tion sn(M, fc) [2J|. The same result can be established 
through one of the infinite product formulas for sn(M, k). 
We arrive at the exact summation of screw dislocations 
for a 7r/2 TGB structure: 



TGB 



(x) = 7z - 



27r 



Im In sn [Ox i^^y, k] 



(18) 



where 9 and -0 are the necessary scale factors, K{k) = 
dx[{l — x^)(l — fc^x^)]~^/^ is an elhptic period, 

iK'{k)^ ^l'^ dx[{l^x^){\-ex'^)]-^/^ ^^K{^/l^) 
is the other elliptic period, and k is the elliptic modulus. 

The ratio of the elliptic periods is r = iK'{k)/K{k). It 
is particularly simple to consider the case that k is pure 
imaginary so that k^ < 0. In this case lmK'{k) — —K{k) 
(see the Appendix), t — Re K'{k)/K{k) + 1 and so we set 
6 = 2K{k)/£d, i/j = ReK'{k)/£t, to achieve the desired pe- 
riodicity. Though we could tune the real part of r, which 
controls the offset between grain boundaries, to generate 
a family of surfaces with different topologies, these iden- 
tities allow for a particularly straightforward analysis of 
the energetics of Schnerk's first surface. Other choices of 
r lead to elliptic moduli in which k is, in general, com- 
plex. 
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The level sets of $tgb, Schnerk's first surface, gener- 
ate the triply-periodic surface shown in Fig.^for £1, — ^d- 
Though Schnerk's first surface is not a niinimizer of the 
smectic free energy, equation it's construction en- 
sures that $TGB is at least a harmonic function. In anal- 
ogy with a single grain boundary which is topologically 
Scherk's first surface, Schnerk's first surface is likewise 
topologically identical to the Schwarz D surface, another 
minimal surface. This can be seen by comparin g th e unit 
cell of Schnerk's first surface (see Fig. 2 of ref. [T3|) with 
that of the Schwarz D surface. Unfortunately, we are not 
aware of a simple transformation, such as a rescaling of 
one or more of the coordinates, that renders Schnerk's 
surface exactly minimal. It may be interesting to note 
that the Schwarz D surface also has a par ametric repre- 
sentation in terms of elliptic functions |25| . 



choices of r (and subsequently k) in terms of an algebraic 
combination of elliptic functions with < 0. Notice, for 
example, that [^^ 



sn(w, ik) 



sn(w\/l + k^, k/y/T- 
dn(u;Vl + fc2, fc/^/^ 



(19) 



allowing us to relate elliptic functions with k^ < to 
those with k^ > 0. Similar relations map cn to cn/dn 
and dn to 1/ dn. Other relations allow us to map elliptic 
functions with modulus k to those with modulus l/Zc, and 
the Landen transforms yield more complicated identities 
between elliptic functions with different moduli. 



IV. PROPERTIES OF SCHNERK'S FIRST 
SURFACE 



B. Schnerk's first surface and the Jacobi elliptic 
functions 



For arbitrary sums of parallel screw dislocations, 
Tjiw) — dx^ — idy^ is a meromorphic function of w whose 
simple poles correspond to screw dislocations. In the case 
of Schnerk's first surface, i] is a doubly periodic function 
of w and we are immediately led to consider generaliza- 
tions of Schnerk's first surface to other defect lattices. 
Fortunately, by Liouville's theorem, a doubly periodic, 
meromorphic function must have residues which add to 
zero^, which in our case means net charge neutrality of 
the screw dislocations. Thus our construction, or any 
generalization of it, must always generate achiral phases 
since the net twist always vanishes. It is only because of 
the special duality at a = 7r/2 that we can construct a 
rotating structure (though the rotation can be thought of 
in either direction). For alternating TGBs at other twist 
angles, adjacent boundaries rotate the layers in oppo- 
site directions. This observation and Liouville's theorem 
provides an alternative derivation of an observation by 
Sethna, found in reference , that constant density con- 
figurations of parallel screw dislocations with the same 
charge are impossible. In fact, a sum of screw disloca- 
tions with the same charge can be performed formally, 
after a suitable regularization, in terms of the Weierstrafi 
elliptic function, C(^) 24']. This function, however, fails 
to be doubly-periodic and r]{w) fails to describe a suit- 
able layer geometry. We will return to the question of 
describing TGB phases with other twist angles in the 
next section. 

We also point out that limiting the discussion to fc^ < Q 
is not too restrictive a condition. By utilizing the vast 
array of elliptic function identities, one can rewrite many 



^ Consider a closed contour around the boundary of one period in 
the complex plane. By periodicity the closed contour integral 
vanishes and thus the sum of the residues must vanish. 



A. Energetics 

To simplify our notation, we define ( — Ox+iipy (not to 
be confused with the Weierstrafi elliptic function), denote 
its complex conjugate as and use Glaisher's notation 
(cs for cn/sn, ds for dn/sn, etc.) for the elliptic functions 
|24| . suppressing the elliptic modulus k. Additionally, 
P'l(C) = PQ(C) since fc^ is real. The compression strain 
for $TGB, Mzz = [1 - (V$)2]/2, is 



l-7'-^(^'+^')|csCdnCf 



87r2 



i/;2) Re [cs^C dn^C] 



(20) 



Because we can choose the periodicity of our structure 
by either altering k (and consequently T)or by altering 
■0 and Q we have some freedom when doing our calcu- 
lations. However, we note that though the symmetries 
are identical as we alter k or (0, V'), the surfaces are 
not. We will discuss this in the following. For now, we 
note that if ^ = B. the compression strain is particu- 
larly simple. For this choice, k is determined through 
i(l - r) = ^<iK'(k)lK{k) = 2lbl£d. 

In the case of a single grain boundary, we set 7 by 
considering y — ±00. Here, the structure is triply pe- 
riodic and we are forced to set |V$| — 1 inside one of 
the periods to set 7. We choose to have the compression 
vanish halfway between the grain boundaries, e.g. along 
y — £iy/2 ov X = £d/4. These lines are also where we 
choose to measure the rotation of the layers, and so this 
is a natural choice. Though we should choose 7 to mini- 
mize the compression energy for a single periodic domain, 
as ib/ td — > 00 these two procedures agree. Our choice of 
k^ < Q is, again, particularly useful in determining 7. On 
the lines C = K{k)/2 + it and C = ^ + jRci^'(fc)/2, for 
i e M, it can be shown that |cs(C, A:) dn(C, /c)|2 = 1 - fc^ 
and thus |V^>| is constant as well. We verify these iden- 
tities in the Appendix. We can, as a result, set 7 along 
these lines and find: 



72 = 1 - b\i - k^)e'^/{2TTf 



(21) 
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Returning to the rotation of the layers, promised in Sec- 
tion we measure the angle from w = £d/'2 — *^6/2 
to u; £d/2 + i£b/2. We find: 



a — 2 sin 



(22) 



Asfc— >0,6'-^7r/€d and these expressions reduce to 
those for a single grain boundary. Requiring a = 7r/2 

I 



sets 7^ = 1/2 and sets £d = b^j2{\ - k'^)K{k)/n. Again, 
for fc^ > 0, or other generic complex values of k, the lines 
of constant tilt angle are no longer straight, making this 
procedure difficult, if not impossible. 



Finally, we make use of the expansion of Jacobi el- 
liptic functions p3l in terms of g = exp [—nK'/K] = 
-exp [-2n£b/ed]: 



In [sn Q = In sin 



TTW 



E 



2'm'KW 

£d 



In 



2q 



1/4 



Vk 



(23) 



r 



to compute the long distance interaction between grain 
boundaries. Recall that in the case of a single grain 
boundary, the nonlinear strain ((SJ leads to power-law in- 
teractions between defects il3| . Here, we are able to ad- 
dress the nature of the interactions between twist-grain 
boundaries. Jacobi's formula H23|) leads to a simple an- 
swer: the interactions are exponential. To see this we 
note that the first term in (|23|l is the lu-dependent part 
of the phase field for a single grain boundary, as in (|12|l . 
Thus, were we to calculate the difference between the 
full TGB structure and a set of non-interacting grain 
boundaries, the correction would come from the infinite 
sum in (|23|1 (note that constant ln(2q^/'*/ Vfc) drops out 
of the energetics). Hence the corrections would be at 
least 0{q) and q — — exp [— 27r^b/^d]: exponentially de- 
caying interactions at £b grows with id fixed. In the 
symmetric case depicted in Fig. 1, though £i, = £d, 
q = — e"^'^ ~= —0.002 is small enough for an expan- 
sion to be reliable. We can compare this result to the 
linear elasticity theory. There, without director modes, 
there is no interaction between screw dislocations. With 
the director modes included, as in screw dislocations 
interact exponentially and it follows that grain bound- 
aries will as well However, the attenuation length in 
this case is the twist penetration length At = ^/K^JB, 
and not the distance between defects, £d- Thus, the ex- 
ponential interactions that we find here are different - 
they arise from elastic strains and not from the optical 
modes of the director. 

There are some complications in calculating the inter- 
action energy which are worth mentioning. To calculate 
the interaction energy we use <I>tgb to evaluate the com- 
pression in one period, £b x id- From this we subtract the 
energy of Ly/£b isolated grain boimdaries, where Ly is the 
dimension of the system in the y direction. However, the 
isolated grain boundaries have tails which extend beyond 
a distance ib- Thus while the intensive interaction energy 



is the difference of two integrals, they are integrals over 
different regions: the first [-£d/2,£d/2] x [-4/2,4/2], 
the second [—id/2, id/2] x [— oo,cx)]. Fortunately, this 
does not spoil our argument that the interaction is a pos- 
itive power of q. In the tail from ^f,/2 to cxd, the integrand 
arising from the single grain boundary falls off exponen- 
tially, also as e~y/^'^ and so the contribution from the tail 
is also proportional to q. We arrive at: 



Ai^c BL, 



2iT£h 



C 



b 
2< 



2 In 



2V2i 



q (24) 



where C is a positive constant of order unity whose pre- 
cise value depends on the choice of cutoff at the dislo- 
cation core, Lz is the z-dimension of the system, and 
an elastic cutoff length ^ is introduced to cutoff a \w\~'^ 
divergence in u^^ near the origin. This cutoff was nec- 
essary in the case of a single grain boundary as well 
and arises from the infinite periodicity at the core of a he- 
lix, i.e. at the core of the defect the spacing between one 
sheet and the next vanishes. Since q = — exp{— 27r4/^(i}i 
this is an attractive, exponential interaction. Stepping 
back and examining the arrangement of defects, the at- 
traction is not at all surprising. The parallel defects in 
adjacent grain boundaries are of the opposite sign and 
would, left to their own devices annihilate. Because of 
this, we conjecture that even when 9 ^, the compres- 
sion energy is attractive. En passant, we note that this 
form of the compression energy allows us to minimize 
over all possible values of the core size. Doing so, we 
find that ^ oc 6, reminiscent of Kleman's split core edge 
defects for large 6 [11 HI . 

We end this section by considering the bending energy. 
For general 9 and ?/', the mean curvature, H — [(1-1- 
^D^yy + (1 + Zy)z^x - 2z^ZyZxy]/{l + + z^ f/"^ is the 
rather foreboding 
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sn2C 



47^2^2 



3/2 



(25) 



r 



where 



i^o = 1 



27r272 

62 



(V'^-6|2)Recs\dn2C 



87r272 



IcsCdnCi' 



and 



47^2^2 



(26) 



(27) 



The mean curvature is finite everywhere and so, when 
calculating the energy of <&tgb : it is not necessary to in- 
troduce a cutoff for the cores. The energy may likewise be 
expanded in powers of q and we again argue for an expo- 
nential interaction. The choice ip = 6 achieved great sim- 
plicity in our analysis of Uzz but the expression for H does 
not suggest a simplifying choice. However, the energy 
calculation does. When subtracting the bending energy 
of the individual grain boundaries, a great simplification 
occurs if the surfaces are minimal, with H ~ 0. Recall- 
ing our discussion of the single grain boundary, if we set 
6 = \/2V' then the fc = limit of H18|l is precisely Scherk's 
first surface, with H = everywhere Thus, the 

interaction energy is positive definite and 0{q^) since H 
is 0{q). We are led to the appealing result that the 
purely repulsive bending energy can balance the attrac- 
tive compression energy and set the preferred value of 

When ^ — 9, a, more natural choice from the point of 
view of the compression energy, the question of the sign 
of the bending energy interaction is much more difficult 
to answer without an involved calculation. We have ver- 
ified numerically, however, that the interaction remains 
repulsive with this choice of scale factors, as well as many 
others, at long distances. In Fig. |2 we plot the bending 
energy for several choices of ratios between tp and 9, span- 
ning from -0 = 9/^/2 to = 9. Only for ip = 9/V2 do 
we find a purely repulsive bending energy for the entire 
range of £b- The cross-over from repulsive to attractive 
bending interaction occurs for £b/id ~ 1, however, and 
we conclude that stable Shnerk phases can exist in cases 
that lb is large enough compared to the dislocation sep- 
aration within a grain boundary. 



B. The triality of Schnerk's first surface 

We have already exploited the duality in the descrip- 
tion of a single twist-grain boundary in constructing 



Schnerk's first surface. Perhaps not surprisingly, Schn- 
erk's first surface not only enjoys the original duality on 
which it was built but, in addition, an added symmetry. 
Put together, Schnerk's first surface enjoys a somewhat 
Kafkaesque triality 29] . The level sets of Schnerk's first 
surface satisfy 



tan(27r7z/6) 



sc[2Kx/ed] _ .sc[iReK'y/£b] 
dii[2Kx/ed] ~ "^dnliReK'y/eb] 



(28) 



We recognize sc(C)/dn(C) as the elliptic generalization of 
tan (, and note that it has a pole at ^ = if of the same 
form as the pole in tan^ at ^. This allows us to view 
the surface as being composed of defects along x, instead 
of z. Near the zeroes in y and z, let Sy = y — 2nib and 
dz = z — mb/(2j) for n,m G Z we have 



sc[2Kx/e4 _ ReK'bSy 
dn[2Kx/£d] 27r74 ^ 

Likewise, near the poles in y and z, let Sy ~ y- 
and 6z = z — {2m + l)(i/4, and 

sc[2Kx/td + K] _ ReK'bSy 
<hi[2Kx/td + K] ^ ~ 2'K^lb Tz 



(29) 

-(2n + l)4 



(30) 



Thus we see that we may view the Schnerk surface as 
made of oppositely charged defects staggered in the yz- 
plane. In fact, this observation recapitulates the duality 
of Scherk's first surface we pointed out earlier. Notice, 
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FIG. 2: Average bending energy as a function of lb /(■<!■ Re- 
sults plotted are for ip — 6/\/2 (diamonds), ip — 6/l.i (stars), 
ip — 6/1.2 (squares), V ~ 0/1.1 (triangles) and ^ = 6 (cir- 
cles); the lines are guides for the eye. For ip = e/^/2, we 
explicitly recover the quadratic dependence on at small q. 
In all cases, the interaction is repulsive for large it /id > 1 
where the small q limit holds. All calculations are for 6 = 2. 
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however, that the rotation (z, x) (a;, —z) coupled with 
b —b does not leave Schnerk's surface invariant (as it 
does for a single TGB). 

Further elliptic gymnastics demonstrates that the un- 
modified Schnerk's first surface may also be viewed as 
a charge neutral arrangement of defects in the xz-plane, 
with cores along the TGB pitch axis. The poles are lo- 
cated at X = {2m + l)£d/2 and z = (2n -I- 1)6/(47) while 
the zeros at x = mid and z = nb/{2j), for m,n E Z. 
Near a zero in x and a pole in z, we have Sx = x — mi^ 
and Sz - z - {2n + 1)6/(47), so that 



. sc[tReK'y/£d] 



Kb 6x 



dnliReK'y/td] TTjid Sz 

Swapping the pole and the zero, we have 6x — x - 
l)£d/2 and Sz = z - nb/{2-/): 



(31) 



{2m- 



T^lid 



5z 



. sc\iReK'y/ld] ^ 

^ dn[iReK'y/ld] (1 - k^) Kb 8x 



(32) 



where k remains the elliptic modulus. Fortuitously, (1 — 
k'^)sc[u]/d-n[u] = d-a[u + i'RcK']/sTL[u + i'RcK'] and so 
becomes 



.sc[iReK'y/ld + iRcK'] _ Kb 5x 
^ dn[iReK'yl£d + iRcK'] ^ ^^n^JI 



(33) 



we note that similar machinations were required in deriv- 
ing 1)30(1 . but we digress. Again we see that we may view 
Schnerk's surface as being built of a lattice of alternat- 
ing defects along the y, or pitch axis. Thus, we observe 
that Schnerk's first surface may be constructed from any 
one of three orthogonal configurations of helicoid-like de- 
fects. We will exploit these bonus dislocations along the 
pitch axis to find alternate topological constructions of 
arbitrary angle TGB phases. 

This discussion suggests utilizing a modified Schnerk 
surface given by the parametric construction |30| 



sgn(6) 



sc[2Kz/e^] sc[2Kx/id 
dn[2Kz/e,] dn[2Kx/ed] " dn\iReK'y/l 



. sc[iReK'y/eh] 
^ ■ (34) 



found by replacing tan(27rz /d) with its elliptic generaliza- 
tion. Here, we have replaced |6|/7 with iz for notational 
simplicity. When = id, as is the case for Schnerk's 
first surface, it is straightforward to see that this surface 
exhibits the correct symmetry with respect to rotating 
the defects by 7r/2 while changing their sign. 

Equation H34() exhibits an additional symmetry which 
we make manifest first by translating the surface along 
the z direction by ^2/2. This gives 



sc[2Kx/id] 
dn[2Kx/id] 



i(l-fc2)sgn(6) 



sc[2A'z/4] sc\iReK'y/ib 



dn[2Kz/iz] dn[iReK'y/ibY 
(35) 

When fc^ = — 1 (fc = i and 2ib — id), elliptic function 
identities yield the equivalent parametric form 

sc[2Kx/id] ^ sc[2Kz/iz] sc[Rei^V/4] , . 
dn[2Kx/id] dn[2Kz/tz]dn[ReK'y/iby ^ ^ 
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FIG. 3: Smectic defects along the pitch axis for a twist angle 
a = 7r/2 — -\/2/32. Positive dislocations are open circles and 
negative dislocations are filled. The circle radius increasing 
with each grain boundary. Two iterations are shown in the 
figure. 



From this it is clear that the translation z z + tz/2 
followed by the rotation (z, y) — > (y, —z) in concert with 
b —b leaves the surface given by equation l|34|) invari- 
ant for iz — ib- Finally, it is also possible to formulate 
this modified surface as a phase field 



^{x,y,z) — sgn(6)Imln 



sc{2Kz/iz,k) 



dn{2Kz/iz,k) 
Imlnsn(^a: -|- iipy, k). 



(37) 



Thus, this new surface can still be thought of as a sum 
of screw dislocations in the xy-plane where we apply a 
nonlinear transformation to the z coordinate rather than 
a simple rescaling by a factor 7. Notice that $ is no 
longer harmonic because of the functional dependence 



V. TOPOLOGICAL CONSTRUCTIONS FOR 
CHIRAL TWIST-GRAIN BOUNDARIES 

The tt/2 TGB phase structure is not chiral, and our 
construction of summing dislocations along the z axis 
(along which the smectic is periodic) necessarily fails 
to produce any chiral structures. We can explore TGB 
phases with twist angles a < 7r/2 by applying an addi- 
tional twist to Schnerk's first surface. From this point of 
view, it is more natural to think of the dual defects along 
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the pitch axis, along y in our conventions, which are com- 
pelled to twist around to follow the surface. If we make 
an analogy between the screw dislocations along the pitch 
axis and a columnar phase of polymers, then the twisting 
of the surface results in the dislocations braiding around 
each other in analogy with the polymers in the columnar 
moire phase [20l |. 

We first briefly review the basic features of the moire 
phase, which arises when polymer chirality competes 
with the columnar order. Much as occurs in the TGB 
smectic phase, chirality can be incorporated into the 
polymer lattice by introducing an ordered arrangement 
of screw dislocations. In reference pO], the authors con- 
sider two cases: either the dislocations form tilt-grain 
boundaries in which entire rows of columns slide past 
each other, or the dislocations form a honeycomb lattice 
which allows groups of polymers to twist around indepen- 
dently. Here however, the lattice we consider is composed 
of two kinds of "polymers" - the positive and negative 
screw dislocations of the smectic. To prevent premature 
aging of the authors and the reader we will refer to the 
original screw dislocations as "smectic dislocations" and 
the screw dislocations in the columnar lattice of the smec- 
tic dislocations as "columnar dislocations" . There will 
not be any point where the distinction will be clear by 
context. 

To simplify our analysis to its most basic level, we 
consider the limit that the grain boundaries are far apart 
so that we may use the k — results for the rotation angle 
of the grain boundaries. The positive smectic defects 
along the pitch axis sit at 



[Znm, ^nni] = [4^, ^<i(m -f i)] 

and the negative at 



(38) 



(39) 



where 4 = b/{2-f) = 6/[2cos(a/2)], 4 = b/[2sm{a/2)] 
(as before), and n,m G Z. For rotation angles a ~ 7r/2, 
— 4 and the smectic dislocations form a pair of in- 
terwoven square lattices with +b smectic dislocations on 
one lattice and —b smectic dislocations on a square lat- 
tice shifted by half a lattice spacing in each direction. 
To form a twist-grain-boundary phase with rotation an- 
gle different from 7r/2, we untwist the structure along 
its pitch axis (i.e. the j/-axis). In doing this we change 
the rotation angle a. At the same time, the smectic dis- 
locations that run along the y-axis get twisted and the 
dislocations must bend to accommodate this mismatch. 
Though the smectic dislocations interact, acquire edge 
components, and suffer any number of other deforma- 
tions, we still expect that when the rotation from one 
slice of the lattice to the next allows some fraction of the 
smectic dislocations to remain straight that there will be 
a local minimum in the free energy. These preferred ro- 
tations will generate a structure analogous to the moire 
phases of chiral polymers. In the case of a = 7r/2, we 
have a perfect agreement of the lattices as the smectic 



dislocations are straight lines along y. For small rota- 
tions away from 7r/2, as shown in Fig. 13 the energy 
grows. Thus for a highly chiral mesogen which favors 
grain boundaries with rotation angles close to 7r/2, this 
local lock-in mechanism could very well force the forma- 
tion of the non-chiral tt/2 grain boundary. Note, how- 
ever, that the nematic director can and will continue to 
rotate with a preferred handedness - it is only the layer 
structure that is achiral. 

At other angles, we must introduce screw dislocations 
into the columnar lattice of smectic dislocations. These 
columnar dislocations will effect the necessary rotations 
of the lattice. In general, if we have a rectangular lat- 
tice with lattice constants iz and id, there will be a 
lock-in angle when we rotate so that the lattice point 
at [{n + i)4, 4/2] rotates onto [{n + i)4, -4/2] , and 



tan(/3„/2) 



(40) 



Note, however, that £d/iz = cot(Q;/2) depends on the 
angle of rotation. Consistency requires that a = Pn and 
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FIG. 4: Smectic dislocations along the pitch axis for a twist 
angle a — tt/S for tilt-grain boundary structure. Positive dis- 
locations are open circles and negative dislocations are filled. 
The circle radius increases with each grain boundary. Shown 
are the dislocations of only two adjacent grain boundaries. 
The solid and dashed lines are columnar dislocations in the 
defect lattice with opposite topological charge. In the next 
pair of adjacent grain boundaries, the columnar dislocations 
rotate by an angle tt/3. The arrows indicate the relative posi- 
tion from the region below the plane of the page to the region 
above the plane of the page. 
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FIG. 5: Smectic dislocations along the pitch axis for a twist 
angle a = 7r/3. Positive dislocations are open circles and 
negative dislocations are filled. The circle radius increases 
with each grain boundary. 



we find 



tan(/3„/2) 



V2r 



r 



(41) 



from wliicli we get the first few moire angles, /3o = 7r/2, 
Pi = 7r/3, and 02 = tan-^l/V^) ~ 48.19°. Remarkably, 
experiments on large angle grain boundary phases |l5ll3l| 
have observed 90°, 60° and 45° grain boundaries, close 
to, if not exactly, the sequence we find here. Note also 
that the precise value of the angles is altered by adjacent 
grain boundaries; the rotation angle we found in l|22() is 



sm 



(I) 



bVl - k^K{k) 



[1 - 4fl - 



(42) 



so that £d 
so that 



(6/[2sin(Q;/2)])(l -4g). This modifies ^ 



tan/?„/2 ; 



V2r 



1 



(43) 



for our structures. The effect is both in the wrong di- 
rection (i.e., making the predicted angles larger) and ex- 
tremely small (e"^'^ « 10"^), so pure geometry is not 
likely to explain the observed twist angles and the (small) 
energetic effects must be included. 

For n = 1, 0„ = 7r/3, which is especially simple because 
the interlocking square lattices becomes a single triangu- 
lar lattice: iz/£d = V^- Though the lattice is invariant 



under rotations by 7r/3, this is only the case when we 
ignore the sign of the smectic dislocations. The handed- 
ness of the rotation and the necessity of columnar dislo- 
cations arises because of the alternating signs. In Fig. 01 
we show the arrangement of columnar dislocations nec- 
essary to rotate the smectic dislocations. Adjacent pairs 
of columnar dislocations have the opposite topological 
charge, so the net charge remains zero from the point of 
view of the undecorated lattice. This arrangement slides 
every other line of smectic screw dislocations relative to 
each other in opposite directions. Though it is not chi- 
ral, chirality appears because, in adjacent columnar grain 
boundaries, the columnar dislocations must rotate along 
with the smectic layers and dislocations. After three of 
these rotations the smectic dislocations return to their 
original arrangement. In Fig. |Slwe show the unit cell of 
the moire phase of the smectic dislocations. 

For n = 2 (a = 2cot~^ VS ~ 48.19°) we consider Figs. 
and As with the usual moire phase of columnar 
liquid crystals, there is a finite fraction of smectic dis- 
locations below the grain boundary which do not match 
up precisely with a smectic dislocation above the grain 
boundary, even ignoring the sign of the smectic disloca- 




FIG. 6; Smectic dislocations along the pitch axis for a twist 
angle a — 2cot~^\/5 (corresponding to n = 2). Positive 
defects are open circles and negative defects are filled. The 
circle radius increases with each grain boundary. Shown are 
the defects of only two adjacent grain boundaries. The solid 
and dashed lines are columnar dislocations in the defect lattice 
with opposite topological charge. In the next pair of adjacent 
grain boundaries, the columnar screw dislocations rotate by 
an angle a. The arrows depict the shifts of the columns from 
below to above the page. 
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FIG. 7: Smectic dislocations along the pitch axis for a twist 
angle a — 2cot~^ y/E (corresponding to n = 2). Positive 
defects are open circles and negative defects are filled. The 
circle radius increases with each grain boundary. Shown are 
the smectic dislocations for four adjacent slabs. 



tion. In Fig. |51we depict a single grain boundary. As 
before, the columnar dislocations in the next slab will ro- 
tate by 2cot~^ y/5. Though there is a coincident subset 
of lattice points in each pair of adjacent grain boundaries, 
this set changes as we move along the pitch axis. In Fig. 
E|we draw one, very complicated, piece of the structure 
composed of just four columnar slabs. As with the tradi- 
tional moire phases, we believe that the allowed rotation 
angles are all irrational fractions of 2tt, but for the cases 
of n = or n = 1. This lack of a repeating pattern 
persists for the higher moire angles giving moire/twist- 
grain-boundary structures of increasing complexity. The 
mind boggles. 

In closing, we conclude that the 7r/3 structure is par- 
ticularly simple from the point of view of the tt/2 struc- 
ture. Having twisted Schnerk's surface to give 7r/3 grain 
boundaries, continuing to change the angle raises the en- 
ergy again from these purely geometric considerations. It 
is interesting to note that some materials in which 7r/2 
TGB phases have been observed also contain regions of 
7r/3 TGB phases |0. Our results suggest that the preva- 
lence of these two rotation angles is due to the apparent 
simplicity of the defect network along the TGB pitch axis 
- other angles presumably raise the energy of the struc- 
ture. 



VI. DISCUSSION 

In summary, we have presented an explicit analytical 
construction for the 7r/2 twist-grain-boundary structure 
by directly summing screw dislocations. This yields a 
phase field from which the energetics of the structure 
can be estimated. Our construction relies crucially on 
a duality: dislocations within one grain boundary can 
be rotated by ninety degrees, as long as the topological 
charge is simultaneously reversed, without changing the 
layers. For a twist angle of precisely 7r/2, the dislocations 
can be made into a parallel, bipartite lattice of positive 
and negative dislocations. 

From this construction emerges another unexpected 
symmetry; we may view the TGB structure as composed 
of a bipartite lattice of defects along the pitch axis. For 
the 7r/2 structure, these defects are straight. Applying 
an additional twist to the 7r/2 structure yields a twist- 
grain-boundary phase with twist angle a < 7r/2. The 
defects along the pitch axis must then twist along with 
the structure, and geometric considerations suggest that 
this raises the energy of the structure. However, we have 
identified a set of twist angles a for which the screw dis- 
locations can be made partially straight, in analogy to 
the moire phase of columnar liquid crystals 20] . The 
structure for a = tt/B is particularly simple, consisting 
of triplets of defects braiding around each other. 

We conclude by noting that Schnerk's first surface is an 
extremum of the linearized smectic free energy. Though 
we evaluated its energy with the nonlinear elastic strain, 
it is not clear how nonlinearities will modify the struc- 
ture of the layers and the interaction between defects. In 
the case of edge dislocations, the defect interactions are 
not strongly modified by the presence of nonlinearities 
(though the layers are) 131 , suggesting that our conclu- 
sions are robust. Finally, we point out that large twist 
angles are likely to result in a decoupling of the director 
from the layer normal. It is also tempting to attempt 
to identify edge dislocations in the smectic with screw 
dislocations in the columnar crystal of defects, and vice 
versa. There is much to do. 



Acknowledgments 

We thank M. Cohen, G. Gibbons, R. Kusner, and T.C. 
Lubensky for interesting discussions. This work was sup- 
ported through NSF Grant DMR05-47230, the Donors of 
the ACS Petroleum Research Fund, and a gift from L.J. 
Bernstein. 



APPENDIX: IDENTITIES FOR fc^ < 

Here we establish some needed identities for the case of 
K pure imagmary, i.e. fc2 e R and fc^ < 0. Though these 
are surely known, we put them here for completeness. 
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The first result we need is Im K'{k) = —K{k). We 
start with the definition of K'{k): 



iK'{k) 



i/fe 



dx 



7(1 



fl - fc2a 



Setting k = iK, for k G R+, we have 



iK'{iK) = 



v/(1-z2)(1 + k2z2) 



(A.l) 



(A.2) 



We choose three branch cuts in the complex plane as 
follows: the first connects —1 to I along the real axis. The 
second starts at i/ n and runs up the imaginary axis, while 
the third starts at —i/n and runs down the imaginary 
axis. We choose a contour from 1 to —i/k which goes 
from 1 to along the real axis and then goes from to 
—i/n along the imaginary axis. We have 



iK'iin) 







dx 



1 ^(1-x2)(1 + k2x2) 

-^/'^ dy 



+ 



(A.3) 



where the overall minus comes from being underneath 
the first cut. Since both integrands are real, it follows 
that ImK'iin) = - dx [{1 - x^){l + K^x^)f'^^ = 
-Kiin). 

Next we show that, again, for k^ E M^, 
|cs(u, k) dn{u, fc)|2 — I — k^ along u — K{k)/2 + it and 
u = t + iKc K' {k)/2. We use the following expressions 
for sn(u), cn(u) and dn(M) in terms of the Jacobi Theta 



functions (using the conventions in |32|) 
sn(w, k) 
cn(M, k) 
dn(u, fc) 





fl, (it 


T ) 


02(O|t) 


ei{v 


t) 


04{O\r) 




t) 


02{0\t) 


0i{v 


t) 


04{0\t) 




t) 


03{0\t) 




t) 



where v — u/[2K{k)]. from which it follows that 

I{u, k) = cs(u, fc) dn(u, fc) 
cn(u, fc) dn(u, fc) 



(A.4) 



(A.5) 



sn(M, fc) 

3|t) e2{v\T)eMr) 



\t) eMT)ei{v\T) 



We have Vl - k^ = 6I|(0|t)/6I|(0|t) ^ so that when fc2 
is real and negative we find 



(A.6) 



Rewriting all the Theta fimctions in terms of 93(v\t), we 
have 



i(u|r)6'3(u|r) 



l-fc-2 



93iv\T)93{v + '^t\t) 



(A.7) 



Fortuitously, 93{v\t) can be expressed in terms of the 
Jacobi triple product: 



03{v\r) = n [1 



^{2'm — l)TTiT~\-2-Kiv 



^(2m — l)7rir — 2x21) 



(A., 



Setting u = ii:(fc)/2 + it {t E M), makes v ^ \ + part. Letting p = e-^^^^^^^'C^^Z-^e^) and g = e"''*/-^(''^) we 
it/[2K{k)\. Since K{k) is real, the real part of v is set. have: 

Moreover, r = 1 + i ^^^(k^^ ^^^'^ ^ prescribed real 



|/(if(fc)/2 + zi,fc)|2 



'9] [1 



ig-i] [l + ip^'^g] [l-ip2mg-i] 



11™=! [1 + ip^'^-^g] [1 - ip2m-l5-l] [1 - ip2rngl^ 



= 1 



(A.9) 



since \a/b\ = \a/b*\. (r - l)/4. Setting h ^ e*'^*/-^^^) we have 

Similarly, when u = t + iReK'{k)/2, v = t/[2K{k)] + 
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I{t + iRcK'{k)/2,k)\^ 


T-rOO 

llm=l 


1 2rn— 

i — p ^ n 




1 2m— f^* 
1 — /9 ^ It 




1 _L T^2m+-^ L 








l-P 


llm=l 
llm=l 


1 

1 -p2m-i^ 




1 +p2m-|/^* 
1 _p2m-|;j* 




1 -p2m+i^ 
1 +p2m+i^ 




1 _p2m-|^* 

1 +p2m-|^* 




llm=l 


1 _p2m-|^ 




1 _p2m+i^* 




1 +p2m-f;j 




1 +p2m-i^* 





1— p2/l*j(l+p 2/1 



1— p 2 h] [1 + h 



= 1 



(A.IO) 



where the final equahty follows from multiplying by 
= 1. Writing C = 9x + iiPy ^ 2K{k)x/td + 
iReK'{k)y/£b, identities (|A.9|I and (|A.10|) show that, 
when 9 = the compression (V$)^ is constant along 



the lines x ~ id/4: and y — plus all their periodic 

translates. These are precisely the lines that run between 
the rows of screw dislocations. 
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